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Abstract 

A new algorithm for computing the multivariate Faa di Bruno's formula is pro- 
vided. We use a symbolic approach based on the classical umbral calculus that turns 
the computation of the multivariate Faa di Bruno's formula into a suitable multinomial 
expansion. We propose a MAPLE procedure whose computational times are faster com- 
pared with the ones existing in the literature. Some illustrative applications are also 
provided. 
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1 Introduction 

The multivariate Faa di Bruno's formula has been recently addressed by the following two 
approaches. Combinatorial methods are used by Costantine and Savits ^Lj, and (only in the 
bivariate case) by Noschese and Ricci |12j . A treatment based on Taylor series is proposed 
by Leipnik and Pearce [10]. We refer to this last paper for a detailed list of references on 
this subject and for a detailed account of its applications. We just mention the paper of 
Savits for statistical applications. A comprehensive survey of the use of univariate and 
multivariate series approximation methods in statistics is given in [S]. 

Computing the multivariate Faa di Bruno's formula by means of a symbolic software can 
be done by recursively applying a chain rule. Despite its conceptual plainess, applications 
of the chain rule become impractical also for small values, because the number of additive 
terms becomes awkward and the derivation of the terms somewhat tiresome. Moreover 
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the output is often untidy, so further manipulations are required to simphfy the result (see 
the example in Appendix 2). So a "compressed" version of the multivariate Faa di Bruno's 
formula becomes more attractive, as the multivariable dimensions or the derivation order 
increase. Here, a "compressed" version of the multivariate Faa di Bruno's formula is given 
by using the umbral methods, introduced and developed in [6l[71[13]. These methods have 
been particularly suited in dealing with topics where the composition of the formal power 
series plays a crucial role, see for instance [H [3] and [5]. Therefore it is quite natural to 
approach the multivariate Faa di Bruno's formula by means of these symbolic tools. The 
result is a new algorithm based on a suitable generalization of the well-known multinomial 
theorem: 



(Xi + X2 H h XnT = ^ 

n 



k\ ko k 



where the indeterminates are replaced by symbolic objects. Suitable choices of these objects 
give rise to an efficient computation of the following compositions: univariate with multi- 
variate, multivariate with univariate, multivariate with the same multivariate, multivariate 
with different multivariates in an arbitrary number of components. 

Finally, the connection between the multivariate Faa di Bruno's formula and the multi- 
nomial theorem allows us to give a closed form for the so-called generalized Bell polynomial 
introduced in [T]. Umbral versions of the multivariate Hermite polynomials are given as a 
special case of these generalized Bell polynomials. 

A MAPLE implementation of all these formulae ends the paper. Comparisons with existing 
algorithms, based on the chain rule, show the improvement of the computational time due 
to the proposed approach. 



2 The umbral syntax 

Umbral methods consist essentially of a symbolic technique to deal with sequences of num- 
bers, indexed by nonnegative integers, where the subscripts are treated as powers. 

More formally an umbral calculus consists of a set A = {a, called alphabet, 
whose elements are named umbrae, and a linear functional E, called evaluation, defined on 
a polynomial ring R\A\ and taking values in R, where -R is a suitable ring. For the purpose 
of this paper, R denotes the real or complex field. The linear functional E is such that 
£^[1] = 1 and 

E[a'p^ ■■■1^] = E[a'] E[fi^] ■ ■ ■ E[-i% (uncorrelation property) (1) 

for any set of distinct umbrae in A and for i,j,k nonnegative integers. A sequence oq = 
1, ai, 02, . . . in i? is umbrally represented by an umbra a when 

^K] = an 

for all nonnegative integers n. The elements {an}n>i are called moments of the umbra a. 
Special umbrae are 
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i) the augmentation umbra e € A, such tha10 E[e^] = (5o,n, for ah nonnegative integers n; 

a) the unity umbra n G yl, such that ii^[n"] = 1, for ah nonnegative integers n; 

Hi) the singleton umbra X ^ ^ such that E[x^] = 5i,n7 for ah integers n > 1. 

Note that a sequence of moments can be umbrally represented by two or more uncorrelated 
umbrae. Indeed, the umbrae a and 7 are said to be similar when 

E[a"] = ^[7"], (in symbols = 7) 

for all nonnegative integers n. An umbral polynomial p is such that p S The support 

of an umbral polynomial is the set of all umbrae which occur. Two umbral polynomials p 
and q are said to be umbrally equivalent if and only if 

E[p] = E[q], (in symbols p c^i q). (2) 

The classical umbral calculus has reached a more advanced level compared with the 
notions that we resume in the next paragraph. We only recall terminology, notation and 
basic definitions strictly necessary to deal with the topic of this paper. We skip the proofs, 
the reader interested in is referred to [H [71 [13] . 

Univariate umbral calculus. The symbol 7'" denotes the product of n uncorrelated 
umbrae similar to 7, that is £'[(7'"')'^] = g^ for all nonnegative integers k and n, where 
gk = E[y^]. The auxiliary umbra 7'"' is called the n-th dot power of 7. The summation of 
n uncorrelated umbrae similar to 7 is denoted by the auxiliary umbra n.7 and called the 
dot product of the integer n and the umbra 7. By using the umbral equivalence ([2|), we are 
able to expand powers of n.7. Indeed, let A = (f^ 2^^,. . .) be an integer partition of length 
vx with multiplicities m(A) = (ri,r2, . . .). Set m(A)! = ri!r2! • • • and A! = (l!)^i (2!)''2 • • • . 
Then powers of n.7 verify the following umbral equivalence: 

where the summation is over all partitions A of the integer z, the symbol {n)y^ denotes the 
lower factorial and 7^ = {li)' ^'^ {12)' '^'^ ' ' ' ■> with 71,72, .. . uncorrelated umbrae similar to 
7. Since E[{n.xy\ = {n)i, equivalence ([3|) can be rewritten as 

(-7)^-E^^^(-xr7A. (4) 

AH ^ ^' 

The main tool of the umbral syntax is summarized in the following construction. By 
equivalence £'[(n.7)*] results to be a polynomial qi{n) = J2\hii^)'^\^>^^['y>^] °^ degree 
i in n, where d\ = z!/m(A)!A!. If the integer n is replaced by an umbra a, the umbral 
polynomial %(«) is such that qi{a) ~ J2xhi(^)'^\d>^'~f>^- denote by a.7 the auxiliary 

^Here Sij denotes the Kronecker's delta. 
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umbra such that (a. 7)* — Qiict), for all nonnegative integers i. The umbra a. 7 is called the 
dot product of a and 7. Since (a.xY — (ck)i, from equivalence we have 

'^^H-E^^^(«-xr7A. (5) 

We repeat this construction, replacing the umbra a in qi{a) with the dot product of a 
and 5, so we construct (a. (5). 7. Since (a. 5). 7 = a.{6.j), parenthesis can be avoided. 

A noteworthy dot product is the auxiliary umbra a.f3.j, where /3 is the so-called Bell 
umbra. The moments of the umbra /3 are the Bell numbers. Indeed, since /3.x = X-f^ = u 
then a. /3.x = Oi.u = a and from equivalence ([5]) we have 

AH "^y^y-^- 

Remark 2.1. Equivalence ([6]) is the umbral version of the univariate Faa di Bruno's for- 
mula. In particular, let {a^}, {gn}-, {^n} denote respectively the moments of 0,7 and a.f3.j. 
If we consider the formal power series 

/(a,t) = 1 + Va„-, /(7,t) = l + y^^n^, /(a./3.7,t) = 1 + V/i^ 

then we have 



n\ ^ — ' n! ^ — ' n! 

n=l n=l n=l 



/(a./3.7,t) = /[a,/(7,t)-l]. (7) 

So £^[(a./3.7)"] is the n-th coefficient of f[a, f{j, t) — 1]. This is why the dot product a./3.7 
is called the composition umbra of a and 7. 



Multivariate umbral calculus. In the univariate classical umbral calculus, the main 
device is to replace a„ with a" via the linear evaluation E. Similarly, in the multivariate 
case, the main device is to replace sequences like {511,22,. where (ii, Z2, . . . , in) G 
is a multi-index, with a product of powers fJ-^i ' ' ' f^n\ where {^i, ^2, . . . , fin} are umbral 
monomials in Note that the supports of the umbral monomials in {/ii, • • • , ^n} 

are not necessarily disjoint. In order to manage a product like ijl^i Ijl^2 ' ' '/^n"' ^ ^ power 
of an umbra, we will use a multi-index notation. Let i € Ng. We set i\ = ii\i2\ ■ ■ ■ in- and 
(^1) Ai2) • • • ) A*n)* = A*iV2^ ■ ■ ■ A^n") where /x = (//i, /i2, • • • , A*n) denotes a n-tuple of umbral 
monomials. We define the i-th power of as follows /i* = fi2 ' ' ' A'rT- 

A sequence {gijieNjj £ -R) with gi = 5ji,i2, and go = 1, is umbrally represented by 
the n-tuple when 

for all i € Nq. The elements {ffijigNff are called multivariate moments of /^. Two n-tuples 
v and /X of umbral monomials are said to be similar, if they represent the same sequence of 
multivariate moments, in symbols v = fi.li for all i,j G Nq, we have E[i/'^fi^] = E[i'^]E[fi^], 
then v and fi are said to be uncorrelated. 

In the univariate case we have defined the auxiliary umbra re. 7 and the composition 
umbra a./S.^. We follow the same steps in the multivariate case. 



4 



Let {fi' , fi" , . . . , fj,'"} be a set of m uncorrelated n-tuples similar to /i. Define the dot 
product of m and as the auxihary umbra m.fj, = fi' + fi" + ■ ■ ■ + fx'" and the m-th dot 
power of fx as the auxihary umbra /x"™' = /x'/x" • • • fi'" . 

In the fohowing we give an expansion of the i-th power of m.fi in terms of partitions of a 
multi-index, see Definition 12.21 To this aim, we introduce the notion of generating function 
of the n-tuple /x. The exponential multivariate formal power series 



gA**^ — gMltl+M2t2H hflntn 



is said to be the generating function of the n-tuple fi, where t = {ti,t2, ■ ■ ■ ,tn) and 
denotes its transpose. Now, assume {gi}ieN^ umbrally represented by the n-tuple If the 
sequence {fi'ijieNjf has exponential multivariate generating function 



k=l i6N; 





by suitably extending the action of E coefficientwise to generating functions ([5]), we have 
-E[e'^*^] = f{fx,t). Henceforth, when no confusion occurs, we refer to f{fx,t) as the gener- 
ating function of the n-tuple /x. 

Proposition 2.1. // {^i, . . . , are uncorrelated umbral monomials, then = jli + ■ ■ ■ + 
where the vectors {p-i, ■ ■ ■ , are uncorrelated and such that /ij = (e, . . . , /ij, . . . , e) is 
obtained from /x by replacing each umbral monomial with e, except the i-th one. 

Proof From ([8]), we have /(/x,t) = E[e'^*'^] = ^[eW*i+/^2*2+-+/^"*"] = n"=i ^[e^'*'] = 
Yll^i f{fii,ti). The result follows by observing that f{fJ,i,ti) = f{fi^,t) for i = 1, 2, • • • , n. 

□ 

In order to compute the coefficients of f{m.fj,,t), we first introduce the notion of com- 
position of a multi-index and then the notion of partition of a multi-index. 

Definition 2.1 (Composition of a multi-index). A composition X of a multi-index i, in 
symbols X \= i, is a matrix X = (Xij) of nonnegative integers and with no zero columns such 
that Ari + Xr2 + ■ ■ ■ + Kk = V for r = 1, 2, . . . , n. 

The number of columns of A is called the length of A and denoted by /(A). 

Definition 2.2 (Partition of a multi- index) . A partition of a multi-index i is a composition 
A, whose columns are in lexicographic order, in symbols X\- i. 

Just as it is for integer partitions, the notation A = (A^^, Ag^, . . .) means that in the 
matrix A there are ri columns equal to Ai, r2 columns equal to A2 and so on, with Ai < 
A2 < • • • . The integer rj is the multiplicity of Aj. We set m(A) = (ri, r2, . . .). 
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Proposition 2.2. Let fi be a n-tuple of umhral monomials. For i € Ng we have 



(9) 



where the sum is over all partitions X of the multi-index i, andfi^ = • • • , 

with . . . uncorrelated n-tuples similar to fi. 

Proof. Due to the uncorrelation property, we have f{m.fj,,t) = [/(/x, t)]"^ and {[/(/i,t) — 
1] + 1}™ = 1 + J2T=i (T) [/(M' *) - I]''- Moreover, we have 

[/(/x,t)-l]'== V ^...^t-i...t-'= 

(7ri,...,7rfc) 

where the sum is over all vectors {tti, . . . , tt^} G Nq \ {0}. If we denote by tt the multi-index 
composition (tti, . . . , tt^) [= i, and we set giTr = Qm ■ ■ ■ 9-Kk ^^'^ (^) = then we have 



i>0 



Replacing ^ in f{m.fi,t) = 1 + (!g^[/(;x,t) - l]^ we have 



E 

!(7r) = fc 



(10) 



f{m.ti,t) = l + Y, 



i>0 



E 

7r|=t 



("^);(7r) 

ttJ /(tt)! • 



fill 



Observe that for each tv \= i, any other composition t \= i, obtained by permuting the 
columns of tt, is such that g„ = Qt. 11 iv = (c^^,C2^, . . .) then the number of distinct 
permutations of tt are Z(7r)!/m(7r)!. So indexing the last summation in (jlip by partitions A 
of i, instead of compositions tt, we have 



/(m./x,t) = l + Yl 



i>0 



(m')/(A) 



.Ahi 



m(A)! A! 



Recalling that EKm.xY^^^] = ("i)«(a) and E[fi^] = gx, equivalence Q follows. 



(12) 



□ 



3 Multivariate Faa di Bruno's formula 



We start by introducing the auxiliary umbra a./3./x. We define 



(a./3./x)^ 



Ahi 



m(A)! A! 



(13) 



where /3 is the Bell umbra and a € A. 
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Theorem 3.1 (Univariate composite Multivariate). Let ^ be a n-tuple of umbral mono- 
mials with generating function f{fx,t). The auxiliary umbra a. 13. fx has generating function 
f{a.l3.ti,t) = f[a,f{ti,t)-l]. 

Proof. The proof follows the same path outlined in the proof of Proposition 12.21 by taking 
backward. Indeed if we set S[a'(^)] = 0;(;!^) and E[nx] = g\, from (HJ]) and equivalence 
(jl3p . the result follows taking into account (jlOp and (jlip and also by observing that 



f{a.(3.fx,t) 



2>0 



9\ 



i\ ^ k\ ^ 

k>l i>0 



7rl— i 



k>l 



TV 



il 



□ 

Theorem 13.11 and equivalence p3]) generalize equation d?]) and equivalence ([6]) respec- 
tively. Indeed the generating function f{a./3.n,t) is the composition of a univariate formal 
power series with a multivariate formal power series and equivalence (jlSp gives its i-th 
coefficient. 

The following theorem characterizes the auxiliary umbra, whose generating function is 
the composition of a multivariate formal power series and a univariate formal power series. 

Theorem 3.2 (Multivariate composite Univariate). Let fi = (^i,...,^„) be a n-tuple of 
umbral monomials with generating function /(/i, t). The auxiliary umbra {^i + • • • + /i„)./3.7 
has generating function 



/[(/il + ■ ■ ■ + //n)./3.7, t] = fifx, (/(7, i) - 1, . . . , /(7, t) - 1)]. 



(14) 



Proof. From Remark |2. 11 we have /[(^i + • • • + fin)-f3-J, t] = E{exp[{ni + • • • + /in)(/(7, t) — 
1)]} = £;{exp(^i[/(7,t) - 1] H h/Un[/(7,t) - 1])}. Hence, the result follows from □ 

The i-th coefficient of /[(/ii + - • ■+ ^n)-f3-l-, t] can be computed by evaluating equivalence 
([6]) after having replaced a with (/ii + • • • + ^„). 



In (inp . if we replace the umbra 7 with the n-tuple u = (z/i, . . . , z/„), we obtain a 
characterization of an auxiliary umbra whose generating function is the composition of a 
multivariate formal power series and a different multivariate formal power series. This 
characterization is given in the following theorem, where the first multivariate formal power 
series is umbrally represented by the n-tuple ^ and the second multivariate formal power 
series is umbrally represented by the n-tuple u. 

Theorem 3.3 (Multivariate composite Multivariate). Let /x = (/Ui, . . . , //„) be a n-tuple 
of umbral monomials with generating function f{n, t) and v = (ui, . . . , Vn) be a n-tuple of 
umbral monomials with generating function f{u, t). The auxiliary umbra {fxi + • • • + iJ.n)-P.i^ 

has generating function 



/[(m + • • • + l^n).f3.U, t] = fi^l, {f{u, t) - 1, . . . , f{v, t) - 1)]. 
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The i-th coefficient of + - • •+ iJ,n)-(3-i^, t] can be computed by evaluating equivalence 
([13]) after having replaced /i with v and the umbra a with {fii + • • • + fin)- 

In Theorem 13.31 we can also choose a m-tuple u = (z^i, . . . , fm), with m 7^ n. To this 
aim, from now on we denote by the vector (ti, ^2; • • • ; tn) of length n. 

Theorem 3.4 (Multivariate composite Multivariate). Let fj, = (fii, . . . , fin) be a n-tuple of 
umhral monomials with generating function f(fj.,t(^n)) '^'^^ = {i^i, ■ ■ ■ ,i^m) be a m,-tuple 
of umhral monomials with generating function f{v',t(^ni))- auxiliary umbra {fii + • • • + 
fin)-f3-i^ has generating function 

/[(/il + • • • + /U„)./3.l^,t(m)] = /[m, (/(l^,t(m)) - 1,- • -Ji^^tim)) " !)]• (15) 

The auxiliary umbra (/ii + - • ■+fin)-f3-i^ is a first umbral expression of the generalized Bell 
polynomials, introduced in [1]. A more general expression can be obtained from equation 
(jl5p . by replacing each occurrence of the umbral monomial u with a m-tuple of umbral 
monomials. In other words, we deal with the composition of a multivariate formal power 
series in n variables and n distinct multivariate formal power series in m variables. 

Theorem 3.5 (Multivariate composite different multivariates). Let fi = {fii, . . . , fin) be a 
n-tuple of umbral monomials with generating function f{p,t(^n)) 0,"^^ {ui,V2, ■ ■ ■ ,i^n} be a 
set of m-tuples of umbral monomials with generating function /(i/j,t(^)), fori = 1,2, • • • ,n. 
The auxiliary umbra fii.fi.Vi + • • • + fin-f^-i^n has generating function 

flfll.p.Ui H h/i„./3.l^„,t(m)] = /[/I, (/(l^l,t{m)) - 1,- • • ,/(l^n,*(m)) " !)]• 

Proof. If we replace the real or complex field R with the polynomial ring R[xi, . . . , x„], then 
the uncorrelation property ([T]) has to be rewritten as [6] 

S[xf • • • • • • ] = • • • E[a']E[P*] • • • , 

for any set of distinct umbrae in A, and for nonnegative integers /c, m, s, t. In i?[xi, . . . , x„] [A\, 
the following property holds 

(xi H h x„)./3.a = (xi./3.a H h x„./3.a). 

In particular we have 

f{xi.p.vi H h x„./3.f„,t) = exp{xi[/(z^i,t) - 1] H h x„[/(i/„,t) - 1]}, 

where {i^i, . . . , f„} are umbral monomials with supports not necessarily disjoint. If we 
replace Vi with the m-tuple Vi, then we still have 

f{xi./3.ui H hXn./3.i^n,*(m)) = exp | ^ Xj [/(i^i, t(^) ) - 1] | . (16) 

Equation (fT6|) still holds if we replace the variables {xi, . . . , x^} with the symbols {/ii, . . . , 

and define the evaluation E on the polynomial ring R[iii, . . . , fin] [A] taking values in , . . . , fin] , 

that is 

f{fll.j3.Ui H h/Un-/3-l^n,*(m)) = exp | ^ /-ij [/(l^j, * (m) ) ~ ' ^^^') 
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Now, if we apply the linear functional E, defined on the polynomial ring . . . , fin] and 

taking values in R, to equation (fT7|) we have 

E[f{fXi.^.Ui H h/X„./3.i^„,t(„))] = flu, (/(i^l,t(m)) - 1, . . . ,/(l^n,t(m)) " 1)], 

by which the result follows. □ 

In [1], the i-th coefficient Bi{xi, . . . , x„) of exp {X^iLi ^Afi^i^ *(m)) ~ M} is called gen- 
eralized Bell polynomial. From equation (|16p . the following corollary is immediately stated. 

Corollary 3.1 (Umbral representation of generalized Bell polynomials). If Bi{xi, . . . ,Xn) 
denotes the generalized Bell polynomial, then we have 

Bi{xi, . . . ,x„) ~ {xi.^.ui H h x„,./3.i^„)\ 

Prom Theorem l3.51 the general expression of the multivariate Faa di Bruno's formula can 
be computed by evaluating the generalized Bell polynomial -Bi(^i, . . . , Hn), as the following 
corollary states. 

Corollary 3.2 (Multivariate Faa di Bruno's formula). If Bi{xi, . . . ,Xn) denotes the gen- 
eralized Bell polynomial, then we have 

Bi{fli, ...,//„) ~ {fJLi.p.Ui H h fln-P-t^nY- (18) 

Note that the umbral polynomial Bi{fii, . . . , fin) is umbrally equivalent to the i-th co- 
efficient of the formal power series exp {X]r=i *(m)) ~ M} i^ CZ])- How to compute 
by means of a symbolic software the "compressed" multivariate Faa di Bruno's formula (jlSp 
is the object of the last section. 

4 Examples and applications. 

Randomized compound Poisson random variable. As observed in [6], the mo- 
ments of a randomized Poisson random variable are umbrally represented by the composition 
umbra a./?. 7. A randomized Poisson random variable is a random sum Sn = ^1 + • • • + A'at 
of independent and identical distributed random variables {-'^j}, where A'^ is a Poisson ran- 
dom variable with random parameter Y. In a.p.j the moments of Y, if they exist, are 
umbrally represented by the umbra a as well as the moments of Aj, if they exist, are um- 
brally represented by the umbra 7. Similarly, the umbra a.(3.fi is the umbral counterpart 
of a multivariate randomized compound Poisson random variable, that is a random sum 
Sn = Xi + ■ ■ ■ + X N of independent and identical distributed random vectors {Xi}. The 
evaluation E, applied to both sides of equivalence (I13p . gives the moments of a randomized 
compound Poisson random variable in terms of moments of {^i} and A^. This result is the 
same stated in Theorem 4.1 of by using a different proof and different methods. Let us 
underline that all auxiliary umbrae, constructed in the previous section, admit probabilistic 
counterpart. 
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Multivariate Laplace transform. In multivariate probability theory, the verifica- 
tion of sign alternation in derivatives of Laplace transform of densities is one of the appli- 
cations for which multivariate Faa di Bruno's formula of all orders is desirable. As it is 
well-known, the Laplace transform Lx{t) of multivariate densities is given by its moment 
generating function, with t replaced by —t, that is Lx{t) = /(a*, — if the moments of the 
random vector X are umbrally represented by the n-tuple /i. Since f{n, —i) = /(— 1./^, t) = 
/(— 1.X./3.//, t), equivalence (|13p can be again used in order to find the derivatives of Laplace 
transform, recalling that E[[ — l.x)^] = {—l)^n\. 

Multivariate cumulants. If x denotes the singleton umbra, multivariate moments 
of x-A* are the so-called multivariate cumulants of /x. Indeed from Proposition 13.11 we have 
f{x-fJ',t) = 1 +log[/(At,t)], by recalling that x-fJ- = X-X-l^-fJ- and f{x-X,t) = 1 +log(t+ 1). 
The equations giving multivariate cumulants in terms of multivariate moments and vice- 
versa are considered in [TT] . See [2] for a recent symbolic treatment of this topic. Here, we 
remark that also equivalence (jlSp allows us to express multivariate cumulants in terms of 
multivariate moments, and vice- versa. Indeed, since x-X-/^-t^ = X-t^^ multivariate cumulants 
in terms of multivariate moments are obtained by applying the evaluation E to both sides 
of equivalence (fT3]l . where the umbra a has to be replaced by the umbra x-X- Analogously, 
multivariate moments in terms of multivariate cumulants are obtained by applying the 
evaluation E to both sides of equivalence (fT3]l . where the umbra a has to be replaced 
by the umbra u and the n-tuple fj, has to be replaced by x-/^) taking into account that 
u.p.{x-tj) = At- 

The multivariate analogous of the well-known semi-invariance property of cumulants is 
stated in Proposition 14.11 First we have to define the disjoint sum /x-i-iv of two n-tuples a* 
and v. As in the univariate case, the auxiliary umbra n+u is such that (fi+uY ~ /i* + 1/*. 
Note that the contribution of mixed products in the disjoint sum is zero. 

Proposition 4.1. // /x and v are uncorrelated n-tuples of umbral monomials, then ^•(a* + 
i") = X-At+X-i^- 

Proof. The result follows by observing that /(/x-j-i/, t) = f{tJ,,t) + f(y,t) — 1. □ 

Thus multivariate cumulants linearize convolutions of uncorrelated n-tuples of umbral 
monomials. By using Propositions 12. l l and 14.1] we are able to state that mixed cumulants 
are zero in the case of uncorrelated components. 

Corollary 4.1. If{f^i}i'^i are uncorrelated umbral monomials, then x-f^ = X'Ai+ ' ' ' +X'An: 

where a* = (/^i, • • • , fJ^n) o-nd the uncorrelated vectors {/ij}"^]^ are obtained from fi by replac- 
ing each component with e, except the i-th one, that is fi^ = {e, . . . , fii, . . . ,e). 

Multivariate Hermite polynomials. The i-th Hermite polynomial Hi{x,T,) is de- 
fined as Hi{x, S) = (— 1)1*1 Di^VC^f^; 0) 5])/(/)(a;; 0, S), where (j){x; 0, S) denotes the multivari- 
ate gaussian density with mean and covariance matrix T, of full rank n. These polyno- 
mials are orthogonal with respect to (^{x; 0, S) if we consider the polynomials Hi(x, T,) = 
S"^), where denotes the inverse matrix of S. The following result is proved 

in [H] 

Hi{x,T.) = E[{xT.-'^ + iYY] Hi{x,^) = E[{x + iZf] (19) 
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where i is the imaginary unit andZ ~ iV(0,S)| and Y ~ A^(0, S ^). By keeping the length 
of the present paper within bounds, we use (fT9]) in order to get the umbral expression of 
multivariate Hermite polynomials. 

Proposition 4.2. If u is a n-tuple of umbral monomials such that f{u,t) = 1 + ^tS^-'^t^ 
and ^ is a n-tuple of umbral monomials such that /(//,t) = 1 + ^tT^t^, then we have 

Hi{x,^)c^{-l.l3.iy + xT,-'^y ^^(a:, S) ~ (-1./3./X + ^e)*. (20) 

Proof Equivalences (|20p follows by observing that from (|19p we have 

°° / 1 \ 

1 + ^ ^ Hi{x,T,)- = expixT.-^t^ - -tT,-H^] , (21) 

k=l iSNjf V / 

°° / 1 \ 

1 + ^ Fi(a;,S)- = expixt^ --tJ:t^] . 

k=l iSNJf ■ ^ / 

\i\=k 

□ 

In ^4], it is proved that Appell polynomials are umbrally represented by the polynomial 
umbra x + a. It is well-known that univariate Hermite polynomials are Appell polynomials. 
Then equivalences (j20p show that also multivariate Hermite polynomials are of Appell type. 
Let us underline that the umbra — 1./3 allows us a simple expression of multivariate Hermite 
polynomials, without the employment of the imaginary unit of equations (119p . Moreover, 
since the moments of Z' ~ N{0,T,) are umbrally represented by the umbra P.f^, then the 
n-th tuple fi = x.(/3./x) umbrally represents the multivariate cumulants of Z. 

The following proposition states that the multivariate Hermite polynomials are special 
generalized Bell polynomials. 

Proposition 4.3. If fi is a n-tuple of umbral polynomials such that f{fi,t) = 1 + ^tS^^t 
and fir^ is a n-tuple of umbral polynomials such that f{fx^,t) = 1 + f{fJ^,t + x) — f{fj,,t), 
thenHi{x,^) ~ (-1)1^1 (-l./3./x^)\ 

Proof. We have 

1 + E E (-i)|^|(-i./3.mj^| = i + E E (-m./xj^^ 

k>l i:\i\=k k>l i:\i\=k 

where —t = (— ti, —t2, ■ ■ ■ , —tn)- The result follows from (12 ip since we have 

1 + E E i?[(-l-/3./^Ji^ = /(-l-/3./^.,-*) = exp('^rS-it^-itS-it^]. 

k>li:\i\=k ^ ^ 

□ 

Finally we remark that computing efficiently the multivariate Hermite polynomials via 
equivalences (j20p or Proposition 14.31 helps in constructing multivariate Edgeworth approxi- 
mation of multivariate density functions, see [9]. 

^ Here Z ~ N{0, E) denotes a multivariate gaussian vector Z with mean and covariance matrix E. 
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5 The UMFB algorithm. 



In this last section, we present a MAPLE algorithm for the computation of the multivariate 
Faa di Bruno's formula by using the umbral equivalence (jlSp . The main steps can be 
summarized as followed: 

i) to the right-hand-side of equivalence (fT8]) . we apply the multivariate version of the 
well-known multinomial theorem: 

(/il./3.l/l H h /in-/3-i^n)* 

^ (^^^_'_^^J(W./3.^l)'^---(^n./3.^nf". (22) 

The procedure mkT finds all the vectors (fci, . . . , such that — 

ii^ Then, it is necessary to expand powers like {iij.p.u The procedure MFB realizes this 
computation by using equivalence (fT3|) , with a replaced by /ij . This procedure makes 
use of the procedure makeTab, available online at 

http : //www.maplesof t . com/applications/view. aspx?SID=33039 

We will add more details on the output of procedure makeTab later on. 

Hi) The procedure joint provides a way to multiply the factors {^i.p.ui)^'^ ■ ■ ■ {^n-P-t^n)^'^ 
in (I22p . previously expanded in ii). 

iv) Finally, occurrences of products like fJ-i^ fJ-'2 ' ' ' /^n" are replaced by gji,j2,...,j„ while occur- 
rences of products like where A = (A^^ , Ag^, . . .), are replaced by (/i^j)^^ i^^x^Y'^ ' ' ' 
where h^^^ denotes the j'-th coefficient of f{ui,t{jf{^). 

All these steps are combined in the procedure UMFB. Some computational results are 
given in Table 1. All tasks have been performed on a PC Pentium(R)4 Intel(R), CPU 2.3 
Ghz, 2.0 GB Ram with MAPLE version 12.0. The procedure dif f of Table 1 is a MAPLE 
routine by which the multivariate Faa di Bruno's formula is computable making use of the 
chain rule. 

Now let us give more details on the quoted procedure makeTab. This procedure gives 
multiset subdivisions. A multiset M is a pair (M, /), where M is a set, called support 
of the multiset, and / is a function f : M ^ Nq. For each fj. e M, the integer f{fj,) is 
called the multiplicity of The notion of multiset subdivision is quite natural and it is 
equivalent to split the multiset into disjoint blocks (submultisets) whose union gives the 
whole multiset. For a formal definition, the reader is referred to pj. As example, for the 
multiset M = {^i, /Ui, /i2}, the subdivisions are 

{{fii,fii,fi2}}; {{/ii,/ii}, {/^2}}; {{^i,/^2, }{w}}; {{w}, {w}, {/^2}}- (23) 

If the input parameter is the vector i, the output of the procedure makeTab gives all the 
subdivisions of a multiset having the vector i as vector of multiplicities, that is f{f^j) = ij for 
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i 


m 


# terms in output 


Time (UMFB) 


Time (diff) 


(6,5) 


2 


14089 


0.7 


1.6 


(7,6) 


2 


60190 


3.2 


29.8 


(7,7) 


2 


123134 


8.1 


75.4 


(5,4) 


3 


20208 


0.7 


2.3 


(6,5) 


3 


122034 


6.3 


62.5 


(5,4) 


4 


86768 


4.5 


26.9 


(5,4) 


5 


288370 


25.9 


130.9 


(4,4,3) 


2 


95138 


6.3 


12.3 


(4,4,4) 


2 


257854 


22.8 


41.1 


(4,3,3) 


3 


313866 


22.5 


54.5 


(4,2,2) 


4 


106912 


6.5 


17.3 



Table 1: Comparison of computational times in seconds. 



J = 1, 2, . . . , n. At a first glance, a multiset subdivision does not seem to have any relation 
with multi-index partitions. But any subdivision of a multiset M, having i as vector of 
multiplicities, corresponds to a suitable partition of the multi-index i. As example, the 
multiset M = {/ii, ni, 112} corresponds to the multi-index (2,1), and the subdivisions in 
(j23|) correspond to the multi-index partitions 




Therefore the procedure makeTab allows us to compute all the partitions of a multi-index 
i. On the other hand, the connection between the combinatorics of multisets and the Faa 
di Bruno's formula was already remarked in [8]. 



6 Appendix 1. 

The UMFB algorithm. 
> 

MFB := procO 
local n,vIndets,E; 
option remember; 
n: =add(args [i] , i=l . .nargs) ; 
if n=0 then return (l);fi; 
vindets : = [seq( alpha [i] ,i=l. .nargs)] ; 
E:=add(f [nops(y[l])]* 
y[2]* 

mul(g[seq (degree (x, vindets [i] ) , 

i=l . .nops (vindets) )] ,x=y [1] ) , 
y=makeTab (args) ) ; 
end: 
> 
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joint := procO 
local pl,p2,M,V; 

V := 'if ' (nargs=l, [[args]], [args] ) ; 
pi := mul (add(y ,y=x) ! , x=ListTools [Transpose] (V) ) ; 
p2 :=mul(x!, x=ListTools ['Flatten'] (V)) ; 
M := max(seq( add(y,y=x) , x=V )); 
expand (pl/p2*mul( eval(MFB(op(V [i] ) ) ,[g=g||i, 
seq(f [j]=f I |i'j,j=l. .M)]), i=l..nargs )); 

end: 
> 

mkT := proc(V,n) 
local vE,L,nV; 
nV : =nops (V) ; 

vE: = [seq (alpha [i] ,i=l. .nV)] ; 

L:=seq( ' if ' (nops (x [1] )<=n,x [1] .NULL) , x=makeTab( op(V) )); 

L:=seq( [seq( [seq(degree(y,z) ,z=vE)] ,y=x) , [0$nV]$(n-nops(x))] ,x=[L] ) ; 

L:=seq( op (combinat [permute] (x)) ,x=[L] ) ; 

end: 

> 

UMFB := procCV, n) 
local S,vE; 
if n=l then 

return (expand (evaKMFB (op (V) ) , [g=gl]))) ; fi; 
vE:=[seq(f | | i=l , i=l . .n)] ; 
S:=add(joint( op(x) ), x= [mkT(V,n)] ) ; 
add( f[ seq(degree(x,f I I i) ,i=l. .n) ]* 
eval(x,vE), x=S ); 

end: 



7 Appendix 2. 

The output of the routine diff of MAPLE, for ■g^^^f{gl{xl,x2),g2{xl,x2)) is 

{£^gl {xl,x2)) (Z)i,i) (/) {gl {xl, x2) , g2 {xl, x2)) £^gl (xl, x2) 
+ i^gl (xl, x2)) (Z)i,2) (/) {gl (xl, x2) , g2 {xl , x2)) £^g2 (xl, x2) 
+ (/) {gl (xl, x2) , g2 {xl , x2)) g^gl (xl, x2) 
+ {£^g2 (xl, x2)) (Z)i,2) (/) {gl {xl, x2) , g2 {xl, x2)) ^gl {xl, x2) 
+ {^g2 (xl, x2)) (1)2,2) (/) {gl {xl, x2) , g2 {xl , x2)) £^g2 {xl, x2) 
+ D2 (/) {gl {xl, x2) , g2 {xl , x2)) Q^g2 {xl, x2). 

The output of the routine UMFB for Q^^^f{gl{xl,x2), g2{xl,x2)) is: 

fhogh,! + f2,o9li,o9io,i + /o,i5^i,i + /o,25^i,o5^o,i + /i,if-^i,of^o,i + h,i9h 
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